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I.  INTRODUCTION 


If  modern  airborne  radiating  systems  are  to  function  properly,  the 
antenna  pattern  must  meet  certain  speci fications .  In  fact,  system 
performance  is  often  very  much  dependent  upon  the  resulting  patterns. 

The  usual  design  procedure  for  an  airborne  antenna  consists  of  allocating 
specific  locations  for  the  antenna  based  primarily  upon  convenience  with 
regard  to  aircraft  structure  specifications.  Measurements  are  then  used 
to  evaluate  the  performance  of  the  antenna  system  in  terms  of  its  desired 
pattern.  This  approach  of  airborne  antenna  design  requires  a  great  deal 
of  engineering  time  and  money.  Thus,  the  need  for  an  analytical  approach 
which  provides  the  antenna  pattern  performance  in  the  presence  of  an  air¬ 
frame  is  quite  apparent. 

One  of  the  approaches  that  has  found  great  success  for  analyzir.q  on- 
aircraft  antennas  is  the  Geometrical  Theory  of  Diffraction  (GTD).  The 
GTD  is  a  high  frequency  asymptotic  solution  in  which  the  scattering 
object  under  consideration  is  large  in  terms  of  wavelength.  In  general, 
a  modern  high-speed  aircraft  is  quite  large  electrically  and  complex  in 
shape.  To  be  able  to  obtain  an  accurate  radiation  pattern,  one  must 
take  the  various  scattering  structures  into  account.  Based  on  its  appli¬ 
cations  to  previous  on-aircraft  antenna  studies  [1-6],  the  GTD  has  proven 
itself  well  suited  to  this  type  of  analytical  study.  Not  only  does  this 
approach  fit  nicely  into  a  ray  optics  format,  but  it  also  provides  a 
means  of  analyzing  the  effects  of  three  dimensional  structures  and  identi¬ 
fying  the  significant  contributions  in  the  resulting  antenna  pattern. 

Previous  GTD  solutions  for  far-field  calculations  have  been  shown  to  be 
accurate  in  predicting  radiation  patterns  when  compared  with  various 
model  measurements.  However,  a  significant  problem  exists  with  attempting 
to  take  far  field  pattern  measurements.  In  order  to  satisfy  the  far-field 
requirement,  one  should  separate  the  transmitter  and  receiver  by  a  minimum 
distance  (i.e.,  ZD^/a,  where  D  is  the  maximum  dimension  of  the  aircraft). 
Using  this  requirement  for  various  scale  models  that  have  been  considered, 
one  would  need  a  range  on  the  order  of  hundreds  to  thousands  of  feet  in 
length.  It  is  obvious,  then,  that  this  requirement  can  not  be  met  for  the 
majority  of  ranges.  The  discrepancy  caused  by  the  near  field  measurement 
is  in  the  definition  of  directions.  That  is,  the  direction  from  the 
antenna  to  the  receiver  is  not  the  same  as  that  from  the  center  of  the 
aircraft  to  the  receiver.  These  directions  are  identical  in  the  far  field. 
Consequently,  real  far  field  measured  data  for  most  airborne  antennas  are 
not  practical.  In  fact  from  a  measurement  point  of  view,  it  is  most  ad¬ 
vantageous  to  measure  patterns  in  the  extreme  near  field  of  the  aircraft 
using  a  small  indoor  range  and  a  scale  model  aircraft.  For  full  scale 
measurements,  it  is  most  convenient  to  make  measurements  of  the  field  in 
the  vicinity  of  the  aircraft  while  the  aircraft  sits  on  the  flight  line. 
Unfortunately,  the  near  field  pattern,  that  is  most  easily  measured,  is 
not  simply  related  to  the  far-ield  pattern.  Therefore,  a  near  field 
analysis  for  the  air-borne  antennas  must  be  carried  out. 
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It  should  be  noted  that  the  GTD  solution  for  the  near-  and  far-field 
analysis  are  simply  related.  Thus  the  concept  of  decomposing  the  aircraft 
structure  into  simpler  components  based  on  previous  analyses  [1-6]  in  the 
far-field  problems  is  still  applicable  in  the  near-field  case.  As  a  first 
step  to  study  the  near  zone  aircraft  problem,  in  this  report  an  analytic 
solution  for  the  near-zone  radiation  from  a  general  elliptic  cylinder  is 
developed.  Section  II  presents  the  high  frequency  asymptotic  analysis  for 
slot  antennas  radiating  from  a  circular  cylinder.  The  results  obtained  in 
Section  II  are  then  generalized,  on  the  basis  of  locality  of  high  frequency 
propagation,  to  elliptic  cylinder  problems  in  Section  III.  Finally,  in 
Section  IV,  numerical  results  calculated  from  the  solution  are  presented 
in  graphical  form  and  compared  with  eigenfunction  solutions.  Excellent 
agreement  is  obtained. 

II.  NEAR  ZONE  H  Ell)S  EXCITED  BY  SHORT  SLOTS  ON 

LARGE  CIRCULAR  CYLINDERS 

Consider  the  radiation  from  a  short  slot  antenna  mounted  on  a  perfectly- 
conducting  circular  cylinder  where  the  surrounding  medium  is  free  space. 

Figure  1  illustrates  the  geometry  of  the  problem,  where  the  antenna  is 
located  at  Q'  and  the  field  point  at  P.  According  to  geometrical  optics 
the  region  exterior  to  the  cylinder  with  the  antenna  at  Q'  is  divided  into 
an  illuminated  region  and  a  shadowed  region  by  a  plane  tangent  to  the 
cylinder  surface  at  Q' .  The  plane  is  referred  to  as  the  shadow  boundary. 

A  portion  of  these  two  regions  adjacent  to  the  shadow  boundary  is  a  transition 
region.  Previous  work  (e.g.,  [7])  based  on  the  Geometrical  Theory  of  Dif¬ 
fraction  (GTD)  formulation  for  the  fields  radiated  from  an  antenna  on  a 
cylinder  employs  three  different  expressions  for  the  fields  in  the  various 
regions.  Namely,  the  geometric-optics  (GO)  solution  is  used  in  the  "deep" 
illuminated  region,  the  creeping-wave  representation  in  terms  of  the  residue 
series  in  the  deep  shadow  region,  and  integral  representation  in  terms  of 
Fock  functions  within  the  transition  region.  One  then  has  to  "blend" 
these  different  representations  for  the  solution  to  obtain  the  total  field 
pattern  over  the  entire  region  exterior  to  the  cylinder.  In  practical  ap¬ 
plications,  this  is  accomplished  by  graphical  extrapolation  of  the  solution 
obtained  using  different  representations  until  they  intersect  smoothly  in  the 
region  where  the  solution  switches  from  one  representation  to  the  other. 

In  this  report  integral  representations  in  terms  of  Fock  functions  are 
employed  throughout  the  whole  space.  The  integral  representation  used 
in  the  entire  illuminated  region  recovers  the  GO  solution  in  the  "deep" 
illuminated  region,  and  the  integral  representation  used  in  the  shadow 
region  reduces  to  the  creeping-wave  (residue)  series  representation  in 
the  deep  shadow  region.  The  following  sections  will  describe  how  the 
high-frequency  near-zone  radiated  field  may  be  calculated  in  the  illumi 
nated  and  shadow  regions.  The  case  for  the  field  point  P  located  in  the 
shadow  region  will  be  investigated  first. 
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A .  Shadow  Region 

Consider  an  infinitesimal  magnetic  current  moment  pm  located  on  the 
surface  of  a  circular  cylinder  at  Q'  as  shown  in  Figure  2.  The  figld 
point  P  is  inside  the  shadow  region.  The  magnetic  current  moment  pm  is 
arbitrarily  oriented  and  may  be  resolved  into  two  orthogonal  components 
in  terms  of  the  two  principle  unit  tangent  vectors  e|  and  e?  of  the  cylinder 
surface  at  Q'  (see  Figure  1).  Thus  an  arbitrarily  oriented  magnetic 
current  moment  can  be  resolved  into  two  orthogonal  components 

Pm  r  *1  Pm  sin  «  +  *2  pm  cos  6  ’  ^ 


where  s  is  the  angle  between  e2  and  pm  as  shown  in  Figure  1.  In  the 
above  expression  p£  and  p$  are  the  equivalent  magnetic  current  moments 
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corresponding  to  a  short  circumferential  and  an  axial  slot  antenna  on 
the  cylinder  surface.  Therefore,  to  represent  the  radiation  due  to  a 

short  slot  of  any  orientation,  it  is  sufficient  to  determine  the  expressions 

for  short,  axial  and  circumferential  slots  mounted  on  a  cylinder.  It  is 

well  known  [8]  that  for  a  short  slot  located  at  Q'  ( a , ^ '  , z ‘  )  on  the 

circular  cylinder  with  radius  a,  the  formal  solution  of  the  fields  external 
to  the  cylinder  at  P(p,<t>,z )  can  be  expressed  as: 

Eaz  =  0  (2) 


ta 


pm  r  -jh(z-z')  1_ 

4,?a  J-  Pp 


(jn)e 


Hn(p'> 


dh 
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where 


9  n  O 

=  k  -  h“. 


The  time  dependence  e^'*  and  the  superscript  for  the  Hankel  function  of 
the  second  kijid  )j . )  are  understood  and  suppressed.  In  the  above 
expressions,  Ea  and  Ec  are  the  fields  due  to  an  axial  and  a  ci rcumferential 
slot  radiating  from  a  circular  cylinder,  respectively.  The  above  formal 
solutions  can  not  be  evaluated  exactly.  However,  for  a  large  radius 
cylinder  and  when  the  field  point  P  is  not  too  close  to  the  cylinder 
surface  (i.e.,  not  inside  the  paraxial  region),  the  above  integrals  can 
be  asymptotically  evaluated.  Following  the  procedure  described  in 
Appendix  I,  the  high  frequency  asymptotic  approximation  (keeping  only  the 
leading  term)  for  the  exact  integrals  can  be  expressed  as: 


Ea  =  0 
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In  the  above  expressions,  g(-)  and  g(-)  are  the  wel 1  - tabul ated  Fock 
functions  [9],  and  the  parameters  a0,  Wq,  siq  and  Sg  are  defined  in 
Appendix  I.  These  vector  components  ot  the  electric  field  can  be  cast 
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into  a  ray  format  in  that  two  orthogonal  field  components  (transverse  to 
the  ray  path)  travel  from  the  source  point  Q',  along  the  geodesic,  to  the 
field  point  P.  This  will  be  described  in  the  following  sections. 

First,  let  us  introduce  several  parameters.  According  to  the 
generalized  Fermat's  principle  proposed  by  Keller  [10],  contribution  of 
F  at  P  is  the  field  on  the  ray  from  Q‘  to  Qi  (which  is  the  geodesic 
on  the  conducting  surface),  then  from  Qi  to  P  which  is  a  straight  line 
tangent  to  the  surface  at  Qi  (see  Figure  2).  The  unit  tangent,  normal, 
and  bi normal  vectors  of  the  ray  are  (tl*  ni  ,  bi=tixn])at  Qi  .  Note  that 
(tl,  n] ,  b] )  form  a  moving  trihedron  along  the  straight  path  from  Qi  to  P. 
The  vector  components  in  terms  of  the  cylindrical  coordinate  system 
(•  ,.,z)  can  be  cast  into  the  vector  components  in  the  ray  coordinate 
system  (t,n,b).  Fquations  (8)-(13),  therefore,  can  be  easily  expressed 
in  the  ray  coordinate  system  as  follows: 

pa 

^  ^  j  2'-  [cos  w1  gUj)  F(«1,s1)n1 


-  cos  g( )  F(^2’s2^n2^  ( ^ 7 ) 
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note  that  n-)  is  the  unit  normal  .vectgr  at,Q-|,  t]  is  the  unit  vector 
directeddfrom  Q]  toward  P,  and  bj  =  t-)  x  n-|  .  Similar  definitions  are 
used  for  (£2*  ^2*  k?)-  ^  's  i n teres t i ng  to  note  that  only  field  com¬ 
ponents  transverse  to  the  ray  exist  in  the  final  expressions.  Before 
proceeding  to  the  next  section,  it  should  be  noted  that  the  final  ex¬ 
pressions,  Equations  (17)- (20 )  are  approximate  solutions  for  large 
radius  cylinders  with  field  point  not  close  to  the  surface.  The  analysis 
for  the  field  point  close  to  the  surface  has  to  be  carried  out  differently. 

B.  Illuminated  Region 

If  the  field  point  P  is  in  the  illuminated  region  as  shown  in  figure  3, 
the  radiation  from  the  source  Q'  on  the  circular  cylinder  propagates  along 
a  straight  line  to  P.  According  to  geometrical  optics,  the  electric  field 
at  P  due  to  a  short  magnetic  current  moment  p  defined  by  Equation  (1)  at 
O'  can  be  expressed  m 


e! 
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where  n'  is  the  unit  normal  vector  at  Q’ ,  t  is  the  vector  directed 
from  Q‘  toward  P,  and  ( 0 1 , ^ ' )  are  the  lofal  ray  coordinate  systems  as 
shown  in  Figure  3.  Note  that  e,  •  t  =  (e,  ■  £)/  sin  9'  and  n'  •  t  = 

(n*  •  s)  sin  0 ' . 

Using  a  heuristic  approach,  Ivanov  [11]  employed  a  parabolic  equation 
method  to  find  the  asymptotic  behavior  of  the  current  in  the  illuminated 
side  of  the  cylinder  due  to  a  plane  wave  illumination.  Ivanov's  solution, 
employing  an  integral  representation  in  terms  of  Fock  functions,  is  valid 
in  the  illuminated  region  and  holds  right  up  to  the  shadow  boundary. 
Furthermore,  his  solution,  when  applied  in  the  "deep"  illuminated  region, 
has  the  same  accuracy  as  the  first  two  terms  of  the  asymptotic  form  of 
the  geometrical  optics  solution.  Based  on  this  observation,  together 
with  the  reciprocity  theorem  and  the  similitude  concept  [12],  Ivanov’s 
result  is  employed  to  modify  the  geometrical  optics. sol ution  given  by 
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Q  ( xs>y$iz$) \  z 


X 


*  / 
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OX - 


—  y 
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Figure  3.  Geometry  relevant  for  the  analysis  in  the  illuminated  region. 


Equations  (21)  to  (23).  The  final  expression  for  the  electric  field  at  P 
inside  the  illuminated  region  due  to  short  slots  on  a  circular  cylinder 
can  he  written  as: 


The  functions  G(-)  and  &(•)  are  related  to  the  Fock  functions  £?(•)  and  o(-) 
by  the  following  equations: 

G(x)  =  e'jx3/3  g(x)  (28) 

and 

G(x)  =  e'jx3/3  g(x)  (29) 

It  should  be  emphasized  that  the  above  Equations  (24)  to  (27)  are  obtained 
via  an  engineering  approach.  This  is  accomplished  by  combining  Ivanov’s 
result  and  the  geometric  optics  solution  for  a  magnetic  current  moment 
located  on  a  large  circular  cylinder.  However,  numerical  results  obtained 
by  using  Equations  (24)  to  (27)  show  excellent  agreement  with  the  eigen¬ 
function  solution  for  the  same  problem  (see  Section  IV),  thereby  con¬ 
firming  their  validity. 
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III.  GTD  GENERALIZATION  FROM  CIRCULAR  CYLINDER  TO 
GENERAL  ELLIPTIC  CYLINDERS 

The  high  frequency  asymptotic  solutions  obtained  for  fields  radiated 
from  slot  antennas  on  a  circular  cylinder  will  now  be  generalized,  on 
the  basis  of  locality  of  the  high  frequency  propagation,  to  the  elliptic 
cylinder  problem.  The  case  for  the  field  point  P  located  in  the  illumi¬ 
nated  region  will  be  investigated  first. 

A.  Illuminated  Region 

Consider  a  slot  antenna  radiating  from  the  surface  of  an  elliptic 
cylinder  with  radius  of  curvature  . 

Based  on  the  local  nature  of  high  frequency  radiation  from  slot  on 
cylinder,  the  results  obtained  in  the  previous  sections  for  the  circular 
cylinder  problem  can  be  directly  applied  to  the  elliptic  cylinder  case. 
Speci f ical ly.  Equations  (24)  through  (27)  can  be  modified  for  the  elliptic 
cylinder  case  by  replacing  the  radius  of  the  circular  cylinder  "a"  by 
the  cross  sectional  radius  of  curvature  p7(Q')  of  the  ellipse  at  the  source 
location  Q1 .  Thus,  the  fields  of  slot  antennas  radiating  from  an  elliptic 
cylinder  can  be  expressed  as: 


lit  e"Jks 
sin  o'  GU]lt)  6  ~ 


COS  0' 


t)  G Um)  e-f --  4. 


(31) 


The  parameters  in  the  above  equations  have  the  same  definition  as  those 
in  Equations  (24)  thorugh  (27).  This  concludes  the  analysis  in  the  il¬ 
luminated  region. 
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B.  Shadow  Region 


In  this  section,  radiation  from  an  elliptic  cylinder  where  the  field 
point  is  located  in  the  shadow  region  is  considered.  Again,  the  solution 
obtained  for  the  circular  cylinder  problem  is  generalized  to  the  elliptic 
cylinder  case  based  on  the  locality  of  high  frequency  propagation. 

One  assumes  that  the  radiation  from  a  cylinder  of  general  elliptic 
cross  section  may  be  described  in  terms  of  rather  highly  attenuated  surface 
waves  originating  in  the  neighborhood  of  the  source  and  traveling  around 
the  cylinder  in  opposite  directions.  For  the  circular  cylinder,  each  of 
these  waves  is  represented  by  a  decaying  amplitude  and  a  phase  that  in¬ 
creases  essentially  linearly  with  distance  along  the  geodesic.  Referring 
to  Equations  (17)  through  (20),  one  notes  that  the  variation  of  the  surface 
wave  is  essentially  described  by  factor  e_Jk*  ■■(  ),  where  <  is  the 
distance  traveled  along  geodesic  path  on  the  surface  and  is  related  to 
7.  The  parameter  <’  given  by  Equations  (15)  and  (16)  is  examined  to 
identify  this  relation.  Equation  (15) 


r. 


1 


(ka  cos  w.j 

"'2' 


(  I  <f  -  4> '  l"a0) 


can  be  written  in  a  different  form: 


1 


f 


1 


(34) 


2 

where  pg  =  a/cos  w-j  can  be  identified  to  be  the  radius  of  curvature 
along  tne  geodisic  path,  and  =  a ( | <f -(f'  | -i0)/cos  wl  is  indeed  the  arc 
length  along  the  geodesic.  It  follows  from  the  work  of  Fock  [13], 

Wetzel  [14],  and  Goodrich  [15]  that  the  procedure  on  an  elliptic  cylinder 
is  to  divide  the  surface  into  small  segments  over  each  of  which  the  radius 
of  curvature  may  be  considered  approximately  constant.  The  variation  of 
the  surface  wave  over  each  segment  is  calculated  as  if  the  segment  were  a 
section  of  a  circular  cylinder  with  a  radius  equal  to  the  local  radius  of 
curvature.  The  total  variation  in  the  surface  wave  is  taken  to  be  the  sum 
of  the  increments  so  calculated.  Thus  the  near-zone  fields,  Equations 
(17)  through  (29),  of  slot  antenna  radiating  from  a  circular  cylinder  can 
now  be  generalized  for  the  elliptic  cylinder  case: 


Ea  * 

*-  _  r\ . 
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j  £  ’I  cos  w!  )^( 


/vv  V/6 
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where  the  integration  is  performed  along  the  geodesic  from  the  source  Q1 
to  the  diffraction  point  .  The  additional  factor 

(V°i>Y/6 

l^Tonj 

must  be  incorporated  in  the  generalization  to  the  general  elliptic 
cylinder  to  preserve  reciprocity  and  to  reduce  uniformly  to  Keller's 
surface  ray  field  in  the  deep  shadow  region  [16]  as  indicated  in  [7]; 
it  is  observed  that  this  factor  also  occurs  in  a  result  given  by 
Logan  and  Yee  [17]  which  is  based  on  a  more  complicated  boundary  layer 
method  of  solution.  This  concludes  the  analysis  of  the  solution  for 
the  near  zone  field  due  to  axial  and  ci rcumferential  slots  radiating 
from  an  elliptic  cylinder. 
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For  an  arbitrarily  oriented  slot  pm  defined  by  Equation  (1),  the 
total  field  can  be  readily  obtained  via  the  superposition  theorem: 

E  -  sin  EC  +  cos  H  Ea  (40) 

n  n  n 

E.  =  sin  r  E^  +  cos  E?.  (41) 

b  b  b 

Note  that  it  is  straightforward  to  transform  the  field  solutions  given  in 
Sections  III  and  IV  into  the  spherical  coordinate  systems  (•■,:)  defined  in 
Figure  3. 

IV.  RICAL  RESULTS 

In  this  section  various  numerical  results  obtained  by  using  the 
equations  described  in  the  previous  sections  are  presented  for  slot 
antennas  radiating  from  a  general  elliptic  cylinder. 

c  *  a 

Consider  a  short  magnetic  current  moment  n|n  -  eipin  sin  f  +  eopm  cos  - 
where  t;  -  45°  radiating  from  a  circular  cylinder.  The  source  is  located 
at  Q'  on  the  cylinder  with  radius  "a"  with  cylindrical  coordinate  system 
(.  '  a,  c ’ - 1 80" ,  and  z'=0).  The  field  point  is  located  at  P  with  (,  ,t,z). 
figures  4  and  5  show  the  patterns  of  Em  and  E;  for  a  field  point  P  located 
at  (.  z  = 1 000 v )  and  f  varies  from  0°  to  360°.  This  is  a  conical  pattern 
with  the  cone  angle  equal  to  45".  The  far-field  conical  pattern  with 
i-45°  is  also  calculated  using  eigenfunction  series  and  plotted  in  Figures 
4  and  5.  It  is  observed  in  Figure  4  that  Eo  calculated  from  the  high  fre¬ 
quency  asymptotic  solutions  at  ,  =1000*  agrees  very  well  with  the  far-/one 
eigenfunction  solution.  The  field  component  E.  as  shown  in  Tigure  5  also 
shows  excellent  agreement  between  the  two  solutions  except  for  340'1- .*  ■  330" . 

The  reason  for  the  disagreement  in  this  small  region  is  that  the  leading  terms 
of  E|  and  E^  given  in  Equations  (17),  (19),  (24)  and  (25)  almost  cancel  with 
each^other  in  that  region.  This  is  demonstrated  in  Figures  6  and  7.  Figure 
6  is  the  pattern  E?  contributed  from  the  circumferential  component  of  the  mag¬ 
netic  current  moment,  and  Figure  7  is  the  pattern  F3  resulted  from  the 
axial  component.  As  described  in  Appendix  I,  for  this  special  situation, 
the  next  higher  order  term  in  the  asymptotic  solutions  is  important  and 
should  be  included.  Figure  8  shows  the  final  result  when  the  next  higher- 
order  term  is  included.  One  observes  that  the  refined  solution  shows 
excellent  agreement  with  the  eigenfunction  result.  Next  some  numerical 
results  for  an  elliptic  cylinder  case  are  presented.  Figures  9  and  10 
show  the  radiation  patterns  Et>  and  E*  due  to  the  same  45"  magnetic  current 
moment  radiating  from  an  elliptic  cylinder.  The  results  seem  to  be 
reasonable,  however,  independent  checks  will  be  performed  later  using 
other  solutions  (e.g.,  an  i ntegral -equat ion  solution).  The  results  just 


shown  confirmed  the  validity  of  the  solution  for  the  far  zone  calculation 
(i.e.,  o=1000A).  However,  the  solutions  are  also  useful  for  near  zone 
calculations.  Figure  11  presents  the  result  E*  in  the  principle  plane 
for  an  axial  slot  radiating  from  a  circular  cylinder  (radius  a-  )  with  the 
field  point  P  located  at  (p=10a,  ,  and  z=0).  The  exact  solution  Equation 

(4)  is  available  but  can  not  be  evaluated  exactly.  However,  numerical 
integration  techniques  can  be  employed  to  obtain  approximate  results.  This 
is  a  tedious  procedure  and  will  be  investigated  in  the  future.  In  order 
to  check  the  near  zone  result,  a  two  dimensional  problem,  nauiely,  an 
infinitely-long  slot  radiating  from  the  circular  cylinder,  is  considered. 

It  is  well  known  that  the  radiation  pattern  for  a  short  axial  slot  in  the 
principal  plane  is  identical  to  that  of  a  magnetic  line  source  mounted  on 
a  cylinder.  The  eigenfunction  result  for  the  two  dimensional  problem  is 
also  plotted  in  Figure  11.  One  can  see  the  excellent  agreement.  For 
comparison,  the  far  zone  eigenfunction  result  is  also  included  in  Figure  11 
It  is  interesting  to  note  that  for  this  zv  ■*.  the  near  zone  pattern  retains 
the  shape  of  the  far  zone  result  but  tiv-  ^  t  to  back  ratio  is  decreased 
in  the  near  zone  case.  Judging  from  ♦  .  m;ults,  it  is  believed  that 

the  solution  obtained  in  Sections  *  , ;■?  a^e  useful  to  calculate  the 
near  zone  field  due  to  short  slot  an:  with  arbitrary  orientation 

radiating  from  a  general  elliptic  c>  'or. 


- -  EIGENFUNCTION  SOLUTION 

X  X  X  X  x  HIGH  FREQUENCY  SOLUTION 

(NO  HIGHER  OROER  TERMS) 

.  0° 

p  =  1000  X  - _ 

Z=I000X  - 

ka  =  10  jr  \ 


-  20 


Figure  4.  Normalized  radiation  pattern  E0  due  to  an  orientated  slot 
antenna  (6=45°)  on  a  circular  cylinder  at  Q'. 
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Figure  5.  Normalized  radiation  pattern  E ^  due  to  an  orientated  slot 
antenna  (b-45°)  on  a  circular  cylinder  at  Q' . 
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Figure  6.  Normalized  radiation  pattern  contributed  from  the  c 
ferential  component  p*-  of  p  in  Figure  5. 
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Figure  7. 


Normalized  radiation  pattern  contrib 
component  p®  of  “pm  in  Figure  5. 
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Figure  10.  Normalized  radiation  pattern  due  to  an  orientated  slot 
antenna  (8=45°)  on  an  elliptic  cylinder. 
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V.  SUMMARY 

Near  zone  radiation  characteristics  for  an  antenna  radiating  from  a 
general  elliptic  cylinder  are  investigated.  Integral  representations  in 
terms  of  Fock  functions  for  the  radiated  fields  are  employed  throughout 
the  entire  space  exterior  but  not  close  to  the  cylinder.  In  the  illumi¬ 
nated  region,  an  engineering  approach,  employing  Ivanov's  results  for  the 
field  produced  by  plane  wave  reflection  at  a  convex  cylinder,  is  used  to 
derive  the  expressions  for  the  fields  due  to  slots  radiating  on  circular 
cylinders.  This  representation  is  valid  in  the  illuminated  region  and 
holds  up  to  the  shadow  boundary.  Furthermore,  this  solution,  when  applied 
in  the  "deep"  illuminated  region,  has  the  same  accuracy  as  the  first  two 
terms  of  the  asymptotic  form  of  the  geometric  optics  solution.  On  the 
other  hand,  a  high  frequency  asymptotic  solution  is  derived  from  the 
rigorous  solution  for  a  slot  antenna  radiating  from  a  circular  cylinder. 

This  solution  is,  again,  an  integral  representation  in  terms  of  the  Fock 
functions.  It  can  be  shown  that  in  the  deep  shadow  region  this  integral 
representation  reduces  to  the  creeping-wave  (residue)  series  representation. 

The  solution  obtained  for  the  circular  cylinder  problem  is  then  gen¬ 
eralized,  on  the  basis  of  locality  of  high  frequency  propagation,  to  the 
general  elliptic  cylinder  case. 

The  solution  described  in  this  report  is  employed  to  calculate  the 
near-zone  field  radiated  from  a  short  slot  mounted  on  a  general  elliptic 
cylinder.  The  numerical  results  thus  obtained  show  excellent  agreement 
with  various  other  solutions  when  the  field  point  is  in  the  far  zone. 

Our  solution  is  also  valid  in  the  near  zone,  however,  it  remains  a  task  to 
generate  the  near  zone  results  by  using  other  solutions  for  comparison. 
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APPENDIX  I.  ASYMPTOTIC  EVALUATION  OF  THE  INTEGRALS 


In  this  appendix  one  of  the  integral  expressions  encountered  in 
Section  II  is  evaluated  asymptotically.  The  procedure  serves  as  a 
recipe,  all  the  expressions  given  in  Equations  (3)  through  (7)  can  be 
evaluated  in  the  same  fashion. 

Consider  now  the  integral 

on 

<  =  ;  e-jh(z'z,)  y  A(n,h)  e'jn(;-r)  dh  (4?) 

n--“- 
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The  first  essential  step  is  to  apply  a  Watson  transformation  to  express 
the  series  in  the  integral  as  a  contour  integral  in  the  manner  [18] 

l  A(n.h)e-jn(‘-*')  -  if  (  A(v.b)  -j»|*-rldv  ,45) 

n=-°°  J  •’c^+c^ 


where  v  is  a  complex  variable,  and  c-|+C2  is  a  contour  which  encloses 
the  poles  of  the  integrand  (i.e.,  the  zeros  of  sin  vu)  as  shown  in 
Figure  12.  Now  the  integral  can  be  expressed  as 


4  - 


1 

2j 


-jh(z-z‘ ) 


A(v,h) 


cl+c2 


gJv-n 

sin  vn 


-Jv  j  i  -<t 


dv  dh  (461 
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Figure  12.  The  integration  contour  in  the  complex  v-plane. 


which  leads  to 


<  =  )'  [T  («}  -  T  (2i-03 

-A  L  in'’  '  m  1 J 
n^O 


where 


T  c , )  ~  r  e'Jh<2-2')  f 


A(-,h)  e 


-jv(  .+2  m) 


d  dh 


■  c. 


i  n  wh  i  c  h  i  -• 


(50) 


(51  ) 


Substituting  Tquation  (43)  into  (51)  and  keeping  only  the  m-0  term, 
one  obtains 


!„(..)  ■  |  e-Jh<2-2,) 


V  ^ 

?  ( ?\ '  e  a  ’  dh 

C1  -a  H  ^  (,a) 


hv 


(52) 


Following  a  standard  procedure  (e.g.  ,  Reference  118]),  one  transforms 
variables  via  the  relation 


fa  V/3 

V  =  ra  +|  2  )  t. 


(53) 


It  follows  from  Wait  [18]  that  in  the  shadow  region,  particularly  in  the 
case  where  |v|  and  Pa  are  both  large  and  of  the  same  order  of  magnitude, 
the  Airy  function  approximation  for  the  Hankel  function  h(^) * (Pa)  is 
given  by 


\2/3  j 

r  /•>  \2/3l  1 

) 

1  •"!'(■,)  J-j 

(54) 


where  wl(t)  is  the  Airy  Integral  [13]  defined  by 


’1 

wUt) 


1 

'/T 


x  e(tx-t3/3)  „„ 


(55) 


The  contour  r ^  is  shown  in  Figure  13. 


30 


Figure  13.  The  integration  contour  in  the  complex  t-plane. 


Returning  to 
6a  %  v;  therefore, 
form 


hJZ)'(6p)  %  ( 


in  which  Bp  cos  a0 
v  ^  pa,  it  follows 
fore 


Equation  (52),  one  notes  that  Bp  are  greater  than 
^vfpp)  ma^  reP^ace^  by  their  Debye  asympotic 


2 _ 

sin  a„ 
0 


-  j  Pp  ( si  na  -a 

sin  aQ  e 


COSa  )  j 

eJ 


IT 

4 


(56) 


2?  2  1 /2 

=  v  so  that  8p  sina  0  =  (r  d  -  v  )  .  Now  since 

that  Bpsinc*0  %  6(p^-a^)l /2  an(j  ao  %  cos-l  (a/p).  There 


sin  «  e 


-jpd 


(57) 


where  d=  („  --a  2)1'2  and  i0  are  shown  in  Figure  1 4 .  Substituting  Equations 
(54)  and  (57)  into  Equation  (52),  the  function  T0 (.* )  can  now  be  written 
as 


(58) 


wj  (t) 


■j) 


d 


si  n 


!o 


-  i .  d  ^  '  'o 
e  '  e 


e  ^ '  1  dv  dh . 


Making  use  of  Fquation  (53),  one  transforms  the  integration  from  -v 
plane  into  t-plane  and  Equation  (58)  is  now  expressed  as 


To(,) 


sin 


1/3 


-  jf  h(z-z  ‘  )  + :- [ d <- a ( >-a  )]> 

0  dh 


(58) 


where 


and 


g(0 


e;ja 
w)  (t) 


dt 


t  e 


-Jtt 


wi(t) 

2 


(60) 


(61) 


(62) 
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Figure  14.  Ray  geometry  in  the  x-y  plane. 


Defining  D  sin  5  =  z-z1  and  D  cos  6  =  d  +  a(.t-^0)  and  transforming  the 
above  integral  TQ(a)  into  w-plane  via  h  =  k  sin  w  [19],  one  obtains 


TQ(a)  %  sin  aQ 


TTsin  w  /.  ,  Y  ka  cos  w  ,  ■ 

-  2  cos  w  1ka  C0SW  +  H  '  2~~~  )  9  U)| 


J 


2j _ 

nkd  cos  w 


- j k D  cos(w-<5)  ,  . 

2  d  '  '  k  cos  w  dw 


(63) 


where  P  is  the  integration  path  in  the  complex  w-plane  shown  in  Figure  15. 

Equation  (63)  can  be  written  in  a  form  suitable  for  an  asymptotic 
evaluation  via  the  method  of  steepest  descent: 


TQ(u)  %  sin  aQ 


f(w)  €kD  q(W)  dw  . 


(64) 


in  which 
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Let  us  now  study  the  geometrical  meaning  of  ws .  Figure  16  illustrates  the 
geometry  of  the  ray  path  from  the  source  Q'  to  the  field  point  P.  It  is 
evident  from  the  developed  view  of  Figure  16(c),  that  the  path  Q'QP 
minimized  the  distance  between  the  source  Q‘  and  the  field  point  P. 

A  part  of  the  ray  path  Q'Q  lies  on  the  cylinder  surface,  where  it  is  a 
geodesic  helix  with  pitch  angle  w$.  The  tangent  line  QP,  then,  completes 
the  ray  path.  Note  that  (see  Figure  16) 


i  cos  w$  =  a(a-ao)  (68) 


s  cos  w$  =  d  .  (69) 

and 

D  =  s  +  i  (70) 


From  the  standard  procedure  (e.g.,  see  [19],  the  asymptotic  approximation 
of  the  integral  Tq ( a )  can  be  expressed  as 


T,W' 


•J 


2tt 

kD[q"(w5Tr 


k0  q(w  ) 

f(ws)e 


(71) 


where 


q(«s)  =  -j 


(72) 


and 


q"(ws)  =  j 


and 


f(ws) 


k  a  sin  w$  gU$) 


■  [i 

.  /sin  w  \  /ka  cos  „  ,1 

*J|(cos^  j(-  r-k)  9  M 


e'jke.  e'jksej  4 


yr 


(73) 


(74) 


35 


(a  )  3D  VIEW 


(c)  DEVELOPED  VIEW 


y 


(b)  2D  VIEW 


Figure  16.  Geometry  for  the  ray  path. 


Substituting  Equations  (72)  through  (74)  into  Equation  (71),  one 
finally  obtains 


T0(a)  * 


|  k2a 


sin  w. 


S(Ss) 


cos  w_ 


\1/3 


g'U5) 


(75) 


e-ikl  e-jks 
Js(s  +0 


where 


/ ka  cos  w  ^  ^ 

f  -  y - J  ( a'a0 )•  (76) 


It  should  be  noted  that  the  second  term  in  Equation  (75)  is  a  higher- 
order  term.  Usually  it  can  be  neglected.  However,  for  some  special 
situations,  (e.g.,  see  Figures  5  and  8  in  Section  IV)  this  higher- 
order  term  is  important  and  should  be  included.  For  consistency,  one 
then  should  also  include  the  next  higher  order  term  for  the  Airy 
function  approximation  of  Hankel  function  given  by  Equation  (56). 
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APPENDIX  II.  LISTING  OF  THE  COMPUTER  PROGRAM  GENERATING 
THE  SOFT  FOCK  FUNCTION  <}S(x) 

r  +j  *3 

6*(x)  =  e  3  a'(x)  x  <  0 

£5 S  ( x )  =  < 


u* ( x )  =  a(x)  x  >  o 


where  a'(x)  is  the  soft  Fock  function  tabulated  by  Logan  [9], 


Kx)  = 


h 


i  xt 

,w 


dt  . 


Note  that  Logan's  table  uses  the  time  dependence  e 


-jwt 


I- 
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ruMPLLX  FUNCTION  IjflIX) 


r  *  **  * 

r  *  *  *  * 
r  *  *  *  * 
r  *  *  *  » 
r  *  ♦  *  * 
C**** 

r.  *  *  *  * 
r  *  *  *  * 
r  *  ♦  ** 
r  **** 


C.b(X)  IS  KELATLO  m  THi  F0CK  FUNCTION  (CURRENT  D  1  S  I  R  I  BU  T  I  UK 
FoMCUON  J  for  I  HE.  SOFT  CASE. 

COP  X<0,0  I . F •  IN  THE  LIT  RFRIONJ 

6R < X ) =CEXP { X**3/3 . *J ) *THE  INTEGRAL  OF 
CEXP»o*X*M/Wl(T  ) /SORT  (PI  ) 

for  x>n.u  i.r.  in  the  shadow  kegion: 

1,H(X)  =  Thl  INTEGRAL  of  CLXP(J+X*T)/kl(T)/SG.KT(PI) 
f'OTE  THAT  THE  Tint  CONVENTION  is  C.t.XPl-J*T> 


for  x>n.u 


PI  PENSION  XF(9>«LL(9)»VFH9).VFD(9).PFL(9»«PF0(9) 

~A1A  II. /t/ 

I'ATA  XF/.0,.b.l.*l.b.?..?.b.3..3,b.4./ 

;  A  T  A  VT  0/. T7bfi. .  AO?9. , ir»7i..f'P6:,,.0?5..('. 191,. 00*3, .00l3b».c/ 

TATA  VFL/. 77SP.1. 377, /.lAut,?.U63. 4, 0/4,5, 01 0,6.004,7. 002, ft.  A  01/ 
n  A  T  A  pFQ/-»,i..  ,-33.  o  7  ,  ,0ft  ob.PH,  71  -  61 . 1  06.  Si  .  14(1 .67  *  170  .  ,  /  00  .  / 
<'ATA  PF  l/-*,0. , -7b. 54, -ft3.p*,-0o.7».-6e.  37, -69,12.-09,48. 
l-dP.67.-f9.7e/ 

PI=3.14it>926t> 

SGN=SIRNC1, .X) 

IFCX.0T.4.)  r-0  TU  /ooo 
TF(X.LT.-4.)  GO  10  3000 
IF(IN.NE.O)  GO  TO  1 
no  ?.  1  =  1.9 
r  l  c I ) =1 . 

DU  2  N=1 . 9 

I  X  =  SG.,*(XF(1  )  -  X  F  ( rj )  ) 

IF (ARS(DX).l T.l.L-6)  0X=1. 

ol«  i)=ci.m*nx 

PONT  I NOE 
I  N  =  1 

OONTINUE 
T  =  l. 

DO  3  N=l«9 
OX=X-XF(N)*SGN 

TF« ABS(DX) .LT.l.E-6)  GO  TO  1000 

t=t*ox 

continue 

VGP=o, 0 
ot,p  =  0  *  0 
nO  4  N  = 1 .9 
VF=VF0<N> 

Pf =PFO<N) 

1F1X.LT. 0.0)  VF  =  VFL(N) 

IF ( X.LT.0.0 )  PF=PFL(N) 

C=CL(N>*(X-XF<N)*SGN) 

r=i ./c 

VGH=VGH+C*VF 

FGR=PGH+c*PF 
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0) 

'4 

ro^TJN'T 

n 

\/(,H  =  Vv»f*  1 

» 

f'bl'rPo  <*  1 

r.b  =  VOH*Cf  XP IlN PL XI  .U,PC^»M/lt'U.  )  ) 

f’t.ruKK 

‘it. 

1  0  l  o 

V  br-  =  Vf  l'1  <  l  ) 

V7 

Pb'liPF1'  <  * » ) 

3' 

IF  (  .  J.C  )  Vf?ti=WFL(.-) 

t)  - 

T  F  (  V  .  LT  .  I.  .  r  )  FC,r(:PFL  (  > 

t: » 

r-or vt-id* rt  xp < r ~,hl>  (  )  i 

b  1 

kt_TUK,., 

b  *•? 

?D0P 

r-d=  (  ,  U  .  .  K  ) 

-  *. 

w^TUK 

’  u  Ci  V 

CHr(.tM.I»v.|P.-(.U,l.  >/(<-.  */**>>) 

-  * 

f  r 

-  >j'i 
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X 


9S(x) 

- ^ - ■> 

majj  phas e  ( deg. ) 


10.000 

20,000 

-69.986 

-9.500 

19.000 

-89.963 

-9.000 

le.ooo 

-89.980 

-C.500 

17. non 

-69,977 

-6.u00 

1  fc  .  0  0  0 

-09.972 

-7.300 

15,000 

-69.966 

-  7  .  u  0  0 

14. OOP 

-89,968 

-fo , 300 

1 3  .  C'  0  0 

-49.948 

-(,.000 

12,000 

-8.9.934 

-b.Df'O 

11.000 

_  #  9|ti 

-  b  .  u  0  0 

1  0  .  r:  0  C 

-89.085 

» *♦ .  o  n  o 

9, con 

-‘'9.643 

-*t .  ui'O 

6.001 

-69.780 

-3.500 

7.002 

-89.670 

-3.000 

6.004 

-B9.48C 

-  2  .  3  0  0 

5.1  10 

-69.120 

-2. UOU 

4.024 

-88,370 

-i .bro 

4.063 

-86.780 

-1  .  C'OO 

2,161 

-f  3.220 

-  ,  bO 0 

1.377 

-75.540 

U.uOU 

.776 

-60,000 

.  bOO 

.363 

-33.370 

1.000 

.167 

.  080 

1.500 

,  066 

35.680 

2  .  Un0 

.  C  25 

71.610 

2  .boo 

.  t  09 

106.510 

4.000 

.003 

140.670 

3 ,  bn  o 

,001 

170,000 

•4  .  UOO 

o.noo 

0.000 

*4.bC0 

o.?oo 

u.ooo 

5.000 

o.noo 

0,000 

b.bOO 

0.  opo 

0.  OOU 

6.000 

0,000 

0.000 

6. 500 

o.noo 

u.ono 

7.0CC 

u.ooo 

U.000 

7.500 

o.noo 

o.noo 

8.UP0 

0.000 

u.ooo 

6.  oPO 

n .  o  o  o 

o.  oon 

9.000 

o.noo 

0,000 

9 .  bO  0 

o.ouo 

0.000 

10.090 

0.000 

o.coo 
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APPENDIX  III.  LISTING  OF  THE  COMPUTER  PRuGRAM  GENERATING  THE  HARD 
FOCK  FUNCTION  gh(x)  and  [yh(x)  J 
f  .  x3 

I  J  3 

G*(x)  =  e  g(x)  x  <  0 

Mh(x)  - 


j  G*(x)  =  d(x) 


x  '  0 


where  q(x )  is  the  hard  Fock  function  tabulated  bv  Logan  [9], 


ixt 


/  ,  lie 

:l(x)  o  dt  • 


Note  that  Logan's  table  uses  the  time  dependence  e 


-jwt 
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1 

2 

3 

u 

S 

t> 

7 

l, 

y 

li* 

n 

IP 

13 

14 

ie> 

i: 

17 

lb 

iy 

2n 

i21 

22 

23 

24 
2b 

2b 

27 

2b 

29 

30 

31 

32 

33 

34 
3b 
3b 
37 
3b 

39 
4  0 

41 

42 

43 

44 
4b 

40 
47 
4B 
49 
b0 


c**»* 

r 

r  **** 
r**** 


r  *  *** 
r  *  *** 

r  **** 


complex  function  g<xj 

r,»X)  IS  RELATED  TO  TH£  FOCK  FUNCTION  l  CURRENT  DISTRIBUTION 
FUNCTION  )  FUR  THc.  HARO  CASE. 

PQR  X<n.U  r.f.  IN  THE  LIT  REGIONS 

G(Xl=C£yP(X**3/3.*J)*TH£  INTEGRAL  OF 
cexp<u*x*t)/wimt)/sqft(pi» 

P  OR  X>0.0  1 . E •  In  the  SHmOOW  nF  G I  UN ! 

G(xj  =  the  integral  of  cexpi  u*x*t  >/h  •  i  t  j/sqkt  c  p  I ) 
A'UTL  THAT  the  time  CONVENTION  IS  CtxPC -uwT  ) 


r  **♦* 

r-l^ENSION  yF(in)*CL(lG)iVFL(10),vFO(ll  ),PFL(10)*PFD(10) 

'ATa  in/u/ 

DATA  VF  0/1  .3994 ,1,232,1 .0591  ,.  73*2,  .4btrl  ,  .31  S3,  .2025  ,.  13, 
1.UH36, . Ub37/ 

DATA  pFQ  /  .  0,4.903,11.36 , 2b.  64, 42.  56,  f  7. 98, 72. 90, 67. 57, 

11 U?. 17, 116.75/ 

DATA  VFL/ 1.3994,1. 552, l.bo2.1 .361,1.948,1 .962,1.994,1  ,99fc. 
21 .9990 .1 ,99°5/ 

DATA  PFL/.f.  , -2*  864, -3. 9, -3. 67i-2.42, -1.45,-. 085«-.  32,-. 33,- 
PI=3. 1415926b 
SijN  =  S  I C-N  t  X  .  .X) 

IF(X.GT,4.)  GO  TO  2000 
TF ( X.LT.-2.5)  GU  TO  3000 
IFUN.NE.O)  GO  TO  1 

no  2  1=1,10 

PL ( I ) =1  . 

DO  ?.  Ns  1,10 

nx  =  1.0MxF(I  )-XF  (N)  ) 

TF ( ARSTOX ) .LT.l.L-6)  L)X  =  1, 

CLCDsCH  J)*OX 
?  continue 

TN=1 

1  CONTINUE 

T  =  l. 

no  3  N=1,10 

rX=X-XF(N),SGN 

IF* ABSTDX) .LT.l.E-6)  GO  TO  1000 
T=T*DX 

3  CONTINUE 

V  G  =  0 . 0 
F&sO.O 
r.u  4  M  =  1 , 1  o 
VF=VFO(N) 

PF=PFQCN) 

TFIX.LT, 0,0)  VFsVFL(N) 

TMX.lt. 0.0)  PF=KFL(N) 

r=CL(N)*(X-XF(N)*RGN) 

r=SGW/f 


22/ 
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31 

vbsvb  +  r.  *vf 

3> 

PG=PG*r*FF 

3  A 

4 

<■  O'lT  Jung. 

3 ‘4 

I 

33 

P  fi  =  H  l ,  ♦  1 

3h 

PsVG*Cf  XP  (  C^PL*  <  .U.Pb*PI/lfcP 

37 

-.t  r  uK-« 

3H 

10^0 

V  o  r  v  F  0  (  INJ  J 

3  .4 

r0rHf  L'(H 

f>  l) 

IF  1  X.t.  1  .1.'  .0  >  I/FSVFLIM 

fcl 

TF«X.LT.r.n»  pG=b  f  l  1  N  I 

b<; 

G  =  \/G*CtXK<i'"PlX(.UiPo*-PI/l£;0 

b  3 

S  4 

?  U  0  0 

-  =  1 .  f'3?s>»CF  x  p  ( rnn  l>  ( 

fc 

Pt.  T  JP  J 

**  r.1 

1  u  o  n 

P(  .0.1.  )/(?.»>»*  3  )-<>./>♦ 

b  / 

f,  =  b-  (  .  Q  .  ^  ♦  7  .  >.  *  »  3 

bf. 

4l  T  JrtU 

by 

run 

...  uoft 

tvr 

•r-^i  i  -  " 


3  *  < ) 
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